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Abstract  Military  recruits  experience  a  high  incidence  of 
febrile  respiratory  illness  (FRI),  leading  to  significant 
morbidity  and  lost  training  time.  Adenoviruses,  group  A 
Streptococcus  pyogenes ,  and  influenza  virus  are  implicated 
in  over  half  of  the  FRI  cases  reported  at  recruit  training 
center  clinics,  while  the  etiology  of  the  remaining  cases  is 
unclear.  In  this  study,  we  explore  the  carriage  rates  and 
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disease  associations  of  adenovirus,  enterovirus,  rhinovirus, 
Streptococcus  pneumoniae ,  Haemophilus  influenzae ,  and 
Neisseria  meningitidis  in  military  recruits  using  high- 
density  resequencing  microarrays.  The  results  showed  that 
rhinoviruses,  adenoviruses,  S.  pneumoniae ,  H.  influenzae , 
and  N.  meningitidis  were  widely  distributed  in  recruits.  Of 
these  five  agents,  only  adenovirus  showed  significant 
correlation  with  illness.  Among  the  samples  tested,  only 
pathogens  associated  with  FRI,  such  as  adenovirus  4  and 
enterovirus  68,  revealed  strong  temporal  and  spatial 
clustering  of  specific  strains,  indicating  that  they  are 
transmitted  primarily  within  sites.  The  results  showed  a 
strong  negative  association  between  adenoviral  FRI  and  the 
presence  of  rhinoviruses  in  recruits,  suggesting  some  form 
of  viral  interference. 


Introduction 

The  unique  susceptibility  of  military  recruits  to  respiratory 
infections,  including  febrile  upper  respiratory  infections, 
pneumonia,  pharyngitis,  and  bronchitis,  has  long  been 
recognized  [1—4].  This  susceptibility  is  assumed  to  stem 
primarily  from  crowding,  and  is  exacerbated  by  stress,  lack 
of  sleep,  and  mixing  of  people  from  diverse  geographic 
locations  [5,  6].  The  primary  recognized  causative  agents  of 
febrile  respiratory  illness  (FRI)  in  recmits  include  adenovirus 
(HAdV),  influenza  virus,  and  group  A  Streptococcus  pyo¬ 
genes  (GAS).  Influenza  and  GAS  are  well  controlled  by 
standard  influenza  vaccination  and  preventive  penicillin 
prophylaxis.  HAdV  is  currently  associated  with  15,000  of 
the  22,000  reported  cases  per  year  of  recmit  FRI  requiring 
hospitalization  or  restriction  to  bed  rest  [7]. 
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Both  GAS  [8]  and  HAdV  [9]  have  been  shown  to  exhibit 
strong  site-specific  strain  maintenance  and  dynamic  strain 
turnover  patterns,  with  specific  strains  dominating  at 
specific  sites  for  periods  of  months  or  years.  For  HAdV, 
the  patterns  of  dominance  at  different  sites  are  essentially 
independent,  such  that  one  serotype  or  genotype  will 
account  for  all  identified  cases  at  one  site,  while  another 
serotype  or  genotype  will  account  for  all  cases  at  another 
site.  This  offers  evidence  that  most  of  the  identified  cases 
result  from  internal  transmission:  infected  recruits  are 
catching  the  pathogen  from  each  other  or  from  the  training 
environment,  rather  than  bringing  the  pathogens  with  them. 

In  addition  to  the  above-mentioned  pathogens,  high  rates 
of  culture-  or  polymerase  chain  reaction  (PCR)-positive 
results  for  potentially  pathogenic  bacteria  and  viruses,  such 
as  Haemophilus  influenzae  and  rhinovirus  (HRV),  among 
otherwise  undiagnosed  FRI  cases  have  led  to  suggestions 
that  these  organisms  might  also  be  transmitting  and 
causing  FRI  among  recruits.  The  role  of  HRV  as  a  potential 
causative  agent  and/or  facilitator  of  infection  remains  an 
area  of  ongoing  research.  Previous  studies  have  indicated 
that  diverse  HRVs  are  present  in  many  populations  [10-12]. 
A  few  studies  also  have  shown  that  HRVs  appear  to  persist 
in  the  respiratory  tract  in  all  ages  without  illness  for  relatively 
long  periods  of  time  [13-15].  Some  viral  infections  of  the 
upper  respiratory  tract  are  known  to  predispose  a  patient  to 
bacterial  superinfection  by  enhancing  bacterial  adherence 
capacity  [16]. 

Certain  strains  of  H.  influenzae ,  Streptococcus  pneumo¬ 
niae ,  and  Neisseria  meningitidis  are  known  to  asymptomat¬ 
ically  colonize  the  nasopharynx  of  substantial  fractions  of  the 
human  population,  and  thereby  are  considered  commensals 
[17-20].  However,  other  strains  of  these  bacterial  species  are 
also  responsible  for  significant  morbidity  and  mortality  due 
to  the  occasional  invasive  infection  [21].  The  carriage  rate 
of  each  species  usually  stays  high  during  childhood  and 
diminishes  in  adulthood  [18]. 

This  study  is  an  exploration  of  the  distribution  patterns 
of  several  bacterial  and  viral  species  among  healthy  recruits 
and  recruits  with  FRI  (both  with  and  without  confirmed 
HAdV  infections).  To  properly  address  carriage  rates  and 
disease  associations  in  a  controlled  manner,  advanced 
technologies  that  can  provide  comprehensive  information 
are  needed.  The  Respiratory  Pathogen  Microarray  (RPM- 
Flu  v.3. 0/3.1  arrays,  Tessarae,  Potomac  Falls,  VA,  USA),  a 
high-density  resequencing  microarray,  can  provide  de  novo 
sequence  data  covering  hundreds  of  thousands  of  base  pairs 
(bp)  from  multiple  targeted  bacteria  and  viruses,  as 
documented  in  previous  studies  [22-27].  RPM-Flu  v.3.0/ 
3.1  arrays  are  designed  to  return  data  from  multiple  genes 
of  86  common  respiratory  organisms  and  high  human 
health  risk  zoonotic  pathogens  likely  to  be  found  in  the 
upper  respiratory  tract,  an  ideal  tool  for  addressing  the 


distribution  and  association  of  multiple  components  of 
the  respiratory  microflora.  The  technology  has  been  success- 
hilly  demonstrated  to  identify  a  broad  range  of  pathogens 
(including  bacteria  and  DNA  and  RNA  viruses)  in  a  single  test 
at  sensitivities  and  specificities  that  are  similar  to  or  improved 
over  those  of  other  technologies  [24,  25] 

In  this  study,  RPM  technology  was  used  to  explore  the 
distribution  of,  and  associations  between,  HAdV,  picoma- 
viruses  (HRV  and  human  enterovirus  [HEV]),  S.  pneumo¬ 
niae ,  N.  meningitidis ,  and  H.  influenzae  among  US  military 
recruits  with  and  without  FRI.  Within  the  limits  of  the 
sample  size,  we  analyzed  the  strain  distributions  of  the 
targeted  organisms.  The  results  show  that  pathogens 
associated  with  recruit  FRI  tend  to  be  characterized  by 
independently  site-specific  strain  circulation,  indicating  that 
they  are  transmitted  primarily  within  sites.  Common 
commensal  organisms  without  clear  disease  associations 
appear  to  be  essentially  randomly  distributed.  Results 
revealed  a  negative  association  between  HRV  and  HAdV, 
as  well  as  between  HRV  and  FRI  (most  of  which  was 
associated  with  HAdV),  suggesting  viral  interference 
between  HRV  and  HAdV.  Results  also  showed  a  positive 
association  among  the  three  most  common  commensal 
bacteria  in  upper  respiratory  tracts,  supporting  observations 
previously  made  in  other  papers  [28]. 

Materials  and  Methods 

Specimen  Collection  and  Sample  Processing 

This  research  has  been  conducted  in  compliance  with  all 
applicable  federal  and  international  regulations  governing 
the  protection  of  human  subjects  in  research.  The  Naval 
Health  Research  Center  (NHRC,  San  Diego,  CA,  USA) 
routinely  collects,  analyzes,  and  archives  throat-swab  speci¬ 
mens  from  consented  recmits  as  part  of  respiratory  infection 
surveillance  at  eight  US  military  training  camps.  If  available, 
approximately  20  recruits  meeting  the  FRI  case  definition 
(oral  temperature  of  >=38°C  and  cough  or  a  sore  throat,  or 
provider-diagnosed  pneumonia)  are  sampled  per  week  at 
clinics  at  each  site,  on  a  walk-in  basis.  Throat  swabs  are 
suspended  in  stabilizing  transport  media  and  archived  frozen 
at  -80°C.  Total  nucleic  acids  were  extracted  from  125 -pi 
throat  swabs  using  the  MasterPure  DNA  and  RNA  Complete 
Purification  Kit  (Epicentre  Biotechnologies,  Madison,  WI, 
USA)  and  dissolved  in  20  pi  of  nuclease-free  water. 

Clinical  specimens  (N=  202)  collected  from  recruits  at 
eight  US  military  training  sites  were  used  for  this  study 
(Fig.  1,  Supplemental  Table  1  and  2).  These  were  divided 
into  two  sets  for  microarray  studies.  In  the  first  set  of  105 
samples,  70  were  collected  during  a  period  of  2  months 
(July  and  August)  in  2007  at  Marine  Corps  Recruit  Depot, 
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Figure  1  Origin  and  clinical 
correlates  of  samples  used  in 
this  study 


San  Diego  (MCRD  San  Diego).  These  70  are  split  into  two 
categories  by  disease  presentation  and  exposure,  35  from 
incoming  recruits  (’’Healthy  Incoming”),  and  35  from 
recruits  in  the  sixth  week  of  training  who  had  not  reported 
for  medical  treatment  (’’Healthy  at  6  Weeks”).  These  were 
compared  with  35  samples  from  recruits  in  training  with 
FRI  symptoms  (’’Typical  FRI”),  collected  between  June  and 
July  of  2007  at  six  different  sites,  including  16  samples 
from  MCRD  San  Diego.  These  105  samples  were  collected 
from  the  same  3 -month  period  and  82%  of  them  were  from 
the  same  site.  Demographic  data  from  this  sample  set  are 
shown  in  Supplemental  Table  1 .  The  second  set,  consisting 
of  97  samples,  was  collected  from  recruits  at  eight  training 
sites  who  reported  with  FRI  over  several  years,  and  were 
chosen  to  maximize  temporal  and  geographic  diversity 
for  strain  distribution  analyses  (Supplemental  Table  2). 
Because  HRV  was  significantly  less  common  than  any  of 
the  bacteria,  we  chose  the  second  sample  set  on  the  basis  of 
previous  positive  HRV  PCR  results.  Hereafter,  these  will  be 
referred  to  as  “HRV+  FRI”  specimens,  though  the  label  is 
not  meant  to  indicate  any  etiological  connection.  This 
sample  set  allowed  us  to  investigate  association  of  HRV 
with  FRI  and  other  viral  and  bacterial  species  and  their 
geographic  and  temporal  distributions  in  a  broad  range  of 
samples. 

Reverse  Transcriptase  and  Polymerase  Chain  Reaction 

HRV  was  screened  in  clinical  samples  using  reverse 
transcriptase  (RT)-PCR  with  the  following  primers  modi¬ 
fied  from  [29]:  F:  5'  GGTCCCATCCCGGAATT  3'  and  R: 
5'  TCCTCCGGCCCCTGAAT  3',  which  amplify  an  approx¬ 
imately  100-bp  sequence  in  the  conserved  5'  untranslated 
region  (5TJTR).  Reactions  were  performed  with  the  OneStep 
RT-PCR  Kit  (QIAGEN  Inc.,  Valencia,  CA,  USA)  with  an 
initial  reverse  transcription  at  50°C  for  30  min,  then 
preliminary  denaturation  at  95°C  for  15  min,  followed  by  40 


cycles  of  94°C  for  15  s,  52°C  for  30  s,  and  72°C  for  30  s. 
Standard  PCR  for  S.  pneumoniae ,  H.  influenzae  [30,  31],  and 
HAdV  [9]  was  performed  as  previously  described. 

Quantitative  real-time  PCR  assays  were  conducted  on 
MyiQ  instruments  (Bio-Rad  Laboratories,  Inc.,  Hercules, 
CA,  USA)  to  determine  the  number  of  viral  genomes  in 
each  sample.  The  findings  for  the  samples  were  compared 
with  those  for  tenfold  serial  dilution  of  prototype  HRV2 
templates  of  known  copy  number  (101  to  106  copies)  by 
using  specific  primers  (HRV-UTR-F1:  TACGGTAACTTTG 
TACGCC  and  HRV-UTR-R1:  ACACGGACACCCAAAG 
TAGT)  and  RT-PCR/PCR  conditions  as  previously  described 
[23].  Real-time  PCR  reactions  were  carried  out  with  iQ 
SYBR  Green  supermix  (Bio-Rad)  in  25- pi  reaction  volumes 
containing  20  mM  Tris-HCl  (pH  8.4),  50  mM  KC1,  3  mM 
MgCl2,  200  qM  each  of  dATP,  dTTP,  dGTP,  dCTP,  SYBR 
Green  I,  10  nM  fluorescein,  200  nM  primers,  0.625  U  iTaq, 
and  4  pi  of  extracted  nucleic  acid. 

Microarray  Processing  and  Analysis 

The  RPM-Flu  v.3. 0/3.1  design  and  specimen  process 
protocols  for  microarray  analysis  have  been  described  in 
detail  in  previous  publications  [22,  27].  Each  sample  was 
subjected  to  one  microarray  assay.  Microarray  processing, 
hybridization,  and  imaging  scanning  were  performed 
according  to  the  manufacturer’s  recommended  protocol 
(Affymetrix,  Inc.,  Santa  Clara,  CA,  USA)  using  a  GeneChip 
Resequencing  Assay  Kit  (Affymetrix)  with  modification  as 
previously  described  [23].  Pathogen  identification  for  the 
RPM-Flu  v.3. 0/3.1  assay  was  performed  using  Computer- 
Implemented  Biological  Sequence  Identifier  Version  2.0 
software  [32],  an  automatic  pathogen  identification  algo¬ 
rithm  based  on  nucleic  acid  sequence  alignment,  which  was 
developed  and  tested  in  detail  in  previous  studies  [24,  27]. 
The  best  match  for  each  consensus  sequence,  as  determined 
through  a  BLAST  search  (word  size  =  7  and  E  value  <10  4) 
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of  the  sequence  against  GenBank  records,  was  regarded  as 
the  most  closely  related  serotype  of  HRV. 

Base  call  rate  was  calculated  as  number  of  base  calls/ 
probe  length  in  each  tiled  sequence.  Heat-map  and  clustering 
dendrograms  were  made  with  dChip  Software  2005  (DNA- 
Chip  Analyzer,  www.dchip.org)  as  previously  described  [27]. 
The  rows  of  the  imported  data  (base  call  rates)  were 
standardized  and  clustered.  Clustering  distance  was  set  at  1 
(correlation  with  average  linkage),  and  gene  ordering  set  by 
cluster  tightness. 

Statistical  analysis  of  base  call  rates  and  titer  for  HAdV 
and  HRV  was  conducted  with  the  Microsoft  Excel  2008 
(Microsoft  Corporation,  Redmond,  WA,  USA)  T  test  func¬ 
tion.  p  values  for  differences  in  observed  vs.  expected 
frequencies  among  pathogen  species  were  calculated  using 
Fisher’s  Exact  binomial  test  [33].  Calculations  were  per¬ 
formed  using  SAS  version  9.1.3  (SAS  Institute  Inc.,  Cary, 
NC,  USA). 


Results 

Microbial  Distribution  in  Typical  FRI  and  Healthy  Samples 

The  first  phase  of  this  study  focuses  on  results  from  35 
Healthy  Incoming  recruits,  35  Healthy  at  6  Weeks  recruits, 
and  35  Typical  FRI  recruits.  The  first  two  categories  in 
combination  are  referred  to  as  “Healthy”.  This  portion  of 
samples  allowed  us  to  generate  disease  association  data  for 
HRV,  S.  pneumoniae ,  H.  influenzae ,  and  N.  meningitidis , 
plus  HAdV  as,  essentially,  a  positive  control  for  etiology. 
As  mild  respiratory  disease  is  ubiquitous  in  recruit  training 
facilities,  some  unreported  subclinical  FRI  may  have  been 
present  in  the  Healthy  at  6  Weeks  population.  The  micro¬ 
array  hybridization  profile  of  these  105  samples  (Fig.  2  and 
Table  1)  showed  three  major  patterns  of  infection  or 
colonization:  HAdV  concentrated  in  Typical  FRI  samples, 
HRV  concentrated  in  Healthy  samples,  and  commensal 
bacteria  (S.  pneumoniae ,  H.  influenzae ,  and  N.  meningitidis) 
common  across  all  samples. 

The  majority  of  Typical  FRI  samples  [71%  (25/35)]  were 
infected  with  HAdV  serotype  4  (HAdV-4),  the  primary 
recognized  cause  of  acute  respiratory  disease  in  basic 
trainees  (Table  1)  [34].  Sequences  generated  from  three 
HAdV-4  tiles  (El A,  fiber,  and  hexon  genes)  were  all  iden¬ 
tified  as  the  currently  circulating  genome  type,  HAdV-4a 
[34],  as  represented  by  GenBank  entry  AY599837.  Further 
analysis  revealed  three  correlated  single  nucleotide  poly¬ 
morphisms  (SNPs),  including  synonymous  nucleotide  sub¬ 
stitutions  in  El  A  and  hexon ,  and  a  nonsynonymous 
substitution  in  fiber  that  clearly  separate  the  identified 
HAdV-4a  viruses  into  two  distinct  genotypes  (genotypes  1 
and  2;  Table  2).  Genotype  2  HAdV-4  viruses  were 


Figure  2  The  microarray  hybridization  profile  of  105  samples 
collected  from  recruits  in  the  categories  Typical  FRI  (FRI),  Healthy 
Incoming  (, H-incoming ),  and  Healthy  at  6  weeks  (H-6  weeks).  Of  35 
Typical  FRI  samples,  16  were  collected  from  MCRD  San  Diego  and 
the  remaining  19  were  collected  at  five  other  training  sites.  All 
Healthy  samples  were  collected  at  MCRD  San  Diego.  HAdV  human 
adenovirus,  HEV  human  enterovirus,  HRV  human  rhino  virus,  H 
Haemophilus  influenzae,  N  Neisseria  meningitidis,  Sp  Streptococcus 
pneumoniae,  Sa  Staphylococcus  aureus 


exclusively  present  in  15  typical  FRI  samples  collected 
from  MCRD  San  Diego.  Genotype  1  HAdV-4a  was 
detected  in  seven  samples  from  Fort  Benning,  Georgia 
and  three  samples  from  Fort  Jackson,  South  Carolina.  Four 
other  samples  generated  hybridization  signals  from  the 
HAdV-7  fiber  gene  tiles.  Resequencing  data  from  these  tiles 
identified  these  as  HAdV- 14,  an  emerging  serotype  of 
HAdV  recently  reported  at  several  recruit  facilities  [9]. 
HAdV- 14  is  not  specifically  tiled  on  the  RPM-Flu  arrays, 
but  can  be  detected  and  identified  from  the  resequencing 
data  yielded  by  tiles  targeted  towards  related  serotypes. 
Besides  HAdV,  two  samples  from  Fort  Benning  were 
positive  for  human  coxsackievirus  A21  (CAV21),  an 
important  cause  of  mild  respiratory  disease  in  military 
personnel  worldwide  [35].  Only  four  HRVs  (one  from  Fort 
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Table  1  Identification  of  respiratory  pathogens  from  trainees  using  RPM-Flu  v.3. 0/3.1 


Category 

Samples 

HAdV-4a 

Type  1 

Type  2 

HAdV- 14 

HRVA 

HRVB 

HEV 

Hi 

Nm 

Sp 

Sa 

FRI 

Othera 

19 

10 

0 

3 

1 

0 

2 

13  (11) 

11 

15  (13) 

1 

MCRD 

16 

0 

15 

1 

2 

1 

0 

15  (15) 

14 

14  (14) 

1 

H-0  wk  MRCD 

19 

0 

1 

0 

6 

0 

0 

14  (14) 

16 

15  (15) 

0 

Pairedb 

16 

0 

0 

0 

5 

0 

0 

13  (11) 

15 

16  (16) 

2 

H-6  wk  MRCD 

19 

0 

1 

0 

6 

2 

0 

18  (18) 

16 

15  (14) 

0 

Pairedb 

16 

0 

3 

0 

8 

2 

0 

13  (13) 

8 

16(15) 

0 

Total 

105 

10 

20 

4 

29 

5 

2 

86  (82) 

80 

90  (87) 

4 

Numbers  in  parentheses  are  samples  identified  positive  by  species-specific  PCR 

HAdV  adenovirus,  HR  V  human  rhinovirus,  HEV  human  enterovirus,  Hi  Haemophilus  influenzae ,  Nm  Neisseria  meningitidis ,  Sp  Streptococcus 
pneumoniae ,  Sa  Staphylococcus  aureus ,  FRI  febrile  respiratory  illness 

a  Samples  collected  from  six  training  sites  other  than  MCRD  (Marine  Corps  Recruit  Depot  San  Diego).  All  other  samples  were  collected  at  MCRD 
b  Paired  “Healthy”  samples  were  collected  from  the  same  16  individuals  at  0  and  6th  week,  respectively 


Benning  and  three  from  MCRD  San  Diego)  were  detected 
in  Typical  FRI  samples  and  all  were  in  samples  co-infected 
with  HAdV-4a.  They  were  identified  as  four  different 
serotypes:  HRV8,  HRV9,  HRV24,  and  HRV84. 

Of  the  Healthy  samples,  41%  (29  of  70)  were  HRV+  (11 
from  Healthy  Incoming  and  18  from  Healthy  at  6  Weeks), 
in  contrast  to  the  11%  HRV+  rate  in  Typical  FRI  samples. 
The  sequence  results  identified  more  than  half  of  them  (18) 
as  most  similar  to  the  HRV1B  serotype.  The  remaining  had 
highest  homologies  to  the  5TJTR  sequences  of  serotypes 
HRV1A  (one),  HRV9  (two),  HRV30  (two),  HRV32  (one), 
HRV45  (one),  HRV100  (one),  and  HRV  serotype  B  (three). 
Among  70  Healthy  samples,  32  were  paired  samples  that 
were  collected  from  16  individuals  at  two  time  points.  In 
this  paired  sample  set,  HRVs  were  identified  in  11  individuals 
either  upon  arriving  at  the  training  site  or  6  weeks  later 
(Fig.  3).  In  particular,  three  incoming  individuals  carrying 
HRV  were  identified  as  having  different  serotypes  6  weeks 
later  (Fig.  3  a).  For  example,  comparison  of  HRV  consensus 
sequences  from  samples  EVR1133T1/T2  and  EVR1146T1/ 
T2  revealed  considerable  divergence  (Fig.  3b),  demonstrat¬ 
ing  turnover  in  colonizing  HRV  serotypes  during  the  first 
6  weeks  of  training. 

About  7%  (5/70)  of  the  Healthy  samples  were  identified 
as  positive  for  genotype  2  HAdV-4a.  Four  of  these  were 
identified  from  recruits  in  the  sixth  week  of  training, 
suggesting  that  these,  like  the  symptomatic  cases  of  HAdV- 


Table  2  Genotyping  of  adenovirus  4a 


E1A-139 

Fiber-436 

Hexon-192 

Genotype  1 

CTG/Leu 

CAT/Arg 

TTT/Phe 

Genotype  2 

TTG/Leu 

CGT/His 

TTC/Phe 

4a,  were  acquired  during  training.  This  rate  of  asymptomatic 
(or  less  symptomatic)  carriage  is  consistent  with  past  reports 
[36].  Opposite  biases  in  the  distribution  patterns  of  HAdV 
and  HRV  in  Typical  FRI  samples  as  compared  to  Healthy 
samples  (Fig.  2)  suggest  that  HAdV  infection,  fever,  or  both 
might  exclude  or  interfere  with  the  presence  of  HRV  (or, 
possibly,  that  HRV  precludes  adenoviral  FRI). 

The  results  from  the  first  105  specimens  also  revealed 
that  all  three  studied  commensal  bacteria  were  common  in 
all  sample  sets  (Fig.  2  and  Table  1).  S.  pneumoniae  was 
detected  in  86%  of  the  samples,  H.  influenzae  in  82%,  and 
N.  meningitidis  in  80%.  Detection  of  high  frequencies  of 
these  commensal  bacteria  among  these  samples  was  further 
confirmed  by  species-specific  PCR  assays  for  H.  influenzae 
and  S.  pneumoniae  (Table  1).  It  was  not  found  that  any  of 
these  bacteria  were  associated  with  either  disease  or  time 
spent  in  the  recruit  facility. 

Distributions  of  HAdV,  HRV,  and  Commensal  Bacteria 
in  HRV+  FRI  samples 

We  extended  the  investigation  of  microbial  flora  to  include 
more  clinical  samples  from  broader  geographic  sites  and 
to  further  understand  the  associations  (or  lack  thereof) 
between  HRV,  HAdV,  and  commensal  bacteria.  To  this  end, 
we  analyzed  97  clinical  samples  collected  from  recruits 
with  FRI  symptoms  at  eight  military  bases  from  2003  to 
2005  (except  for  two,  one  each  from  2000  and  2006),  all  of 
which  were  previously  identified  as  HRV+  by  PCR  [29] 
(Supplemental  Table  2).  We  also  analyzed  five  cultured 
HAdV- 14  samples  and  one  blank  used  as  blind  negative 
controls  for  HRV. 

Table  3  and  supplemental  Figure  1  show  that  89%  (86/97) 
of  the  HRV+  FRI  samples  were  HRV+  and/or  HEV+  by 


4?)  Springer 


628 


Z.  Wang  et  al. 


Figure  3  a  Microarray  hybrid¬ 
ization  profile  of  32  paired 
samples,  collected  from  16 
Healthy  individuals  on  arrival 
and  again  6  weeks  after  arriving 
at  the  MCRD  San  Diego  train¬ 
ing  site.  In  each  paired  sample, 
the  upper  row  (Tl)  represents 
Healthy  Incoming  and  the  lower 
row  (T2)  represents  Healthy  at 
6  Weeks.  Recognized  strains  of 
human  rhino  virus  ( HRV)  most 
closely  matching  the  clinical 
strains  are  designated  in 
parenthesis,  b  Two  pairwise 
comparisons  of  HRV  consensus 
sequences  from  the  same  indi¬ 
vidual  collected  at  different  time 
points.  Nucleotide  discrepancies 
are  framed  in  boxes 


HAdV  HEV 


HRV 


H  NSp 


ti — if- 


f ; 


1102T1 
1102T2 
1105T1 
1105T2 
1107T1 
1107T2 
1109X1 
1109T2 
1117T1 
1117T2 
1125T1 
1125T2 
1127X1 
1127T2 
1133X1 
1133X2 
1134T1 
1134X2 
1141T1 
1141T2 
1143X1 
1143T2 
1146T1 
1146X2 
1147T1 
1147T2 
1148X1 
1  148X2 
1150T1 
1150T2 
11S5T1 
1155T2 


(HR  VI 00) 


(HRV9) 
(HR  VI  A) 
(HR  VI B) 


[HR  VI B) 
(HRV32) 
(HRVB) 
(HR  VI B) 
(HRV30) 

(HR  VI B) 

(HRV1B) 

(HR  VI B) 
(HR  VI B) 


-4i  -5 *  -SJ  -5.1  -t.T  -JJ  -U  -I 


B 


RPM-Flu  analysis,  and  51%  (49/97)  were  HAdV-4+.  The 
colonization  rates  of  commensal  bacteria  in  this  sample  set 
were  H.  influenzae  (88%,  85/97),  S.  pneumoniae  (81%, 
79/97),  and  N.  meningitidis  (76%,  74/97)  (Table  3),  similar 
to  the  rates  seen  in  the  initial  sample  set.  In  addition,  a  few 
other  common  respiratory  viruses  and  bacteria  (influenza 
A/H3N2,  parainfluenza  3,  Pseudomonas  aeruginosa ,  Klebsiella 
pneumoniae ,  Moraxella  catarrhalis ,  Mycoplasma  pneumoniae , 
and  S.  pyogenes)  were  occasionally  detected  (Table  3).  Overall, 
the  majority  of  samples  were  colonized  with  three  to  six  of  the 


respiratory  viruses  and/or  bacteria  detectable  with  the  RPM-Flu 
v.3. 0/3.1  assay. 

Adenovirus 

Detailed  analysis  of  HAdV-4  resequencing  information 
generated  on  the  RPM-Flu  v.3. 0/3.1  arrays  from  the  49 
HAdV-4+  HRV+  FRI  samples  resulted  in  identification  of 
two  different  strains:  HAdV-4p4  (an  HAdV-4  prototype  or 
vaccine  strain,  only  recognized  as  circulating  at  the  Coast 
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Table  3  Summary  of  respiratory  pathogen  identification 

in  eight  training  sites  with  RPM-Flu  v.3. 0/3.1  microarrays 

Location 

Sample 

Date 

HAdV 

HRV 

HEV 

Hi 

Spn 

Nm 

Mp 

Pa 

Kp 

Spy 

Me 

Flu 

PIV3 

4 

14 

2 

A 

B 

MCRD-SD 

37 

2003-2005 

22 

0 

1 

27 

3 

7 

34 

30 

31 

2 

5 

3 

2 

1 

MCRD-PI 

19 

2003-2006 

12 

2 

1 

15 

2 

0 

18 

14 

14 

1 

3 

1 

4 

1 

Lackland 

AFB 

14 

2000-2005 

3 

0 

10 

0 

0 

10 

14 

9 

1 

4 

1 

1 

Ft.  Jackson 

10 

2003-2005 

4 

0 

8 

0 

0 

8 

8 

8 

1 

2 

Ft.  Leonard 

7 

2004-2005 

4 

0 

4 

0 

0 

6 

4 

5 

1 

Wood 

Ft.  Benning 

6 

2003-2005 

1 

1 

5 

0 

1 

5 

5 

4 

1 

1 

CGTC 

3 

2003-2005 

3 

0 

3 

0 

0 

3 

3 

2 

1 

1 

NSTC 

1 

2005 

0 

0 

1 

0 

0 

1 

1 

1 

0 

1 

0 

0 

Total 

97 

2003-2006 

49 

3 

2 

73 

5 

8 

85 

79 

74 

4 

12 

9 

6 

6 

1 

1 

HAdV  adenovirus,  HR  V  human  rhino  virus,  HEV  human  enterovirus,  Hi  Haemophilus  influenzae ,  Spn  Streptococcus  pneumoniae ,  Nm  Neisseria 
meningitidis ,  Mp  Mycoplasma  pneumoniae ,  Pa  Pseudomonas  aeruginosa ,  Kp  Klebsiella  pneumoniae ,  Spy  Streptococcus  pyogenes ,  Me  Moraxella 
catarrhalis,  Flu  influenza  virus,  F7F  parainfluenza  virus,  MCRD  Marine  Corps  Recruit  Depot,  SD  San  Diego,  P/  Parris  Island,  CGTC  Coast 
Guard  Training  Center,  Cape  May,  NJ,  NSTC  Naval  Service  Training  Command,  Great  Lakes,  IL 


Guard  Training  Center  Cape  May,  New  Jersey,  in  the  last 
10  years  [34])  and  HAdV-4a  (common  at  all  other  sites 
[34]).  The  two  samples  with  HAdV-4p4  were,  in  fact, 
collected  at  Cape  May,  while  the  remaining  47  samples 
were  collected  at  other  sites.  Further  analysis  of  the  HAdV- 
4a  strains  based  on  SNPs  in  fiber ,  hexon ,  and  El  A  genes 
revealed  two  genotypes  matching  those  in  the  Typical  FRI 
sample  set:  genotype  2  HAdV-4a  was  only  present  at  the 
MCRD  San  Diego  site  while  genotype  1  HAdV-4a  was 
circulating  at  the  remaining  training  sites.  Hence,  all  three 
detected  HAdV-4  variants  were  geographically  and  tempo¬ 
rally  segregated,  with  strong  agreement  between  the  sample 
sets  and  between  the  data  generated  here  using  the  RPM- 
Flu  v.3. 0/3.1  assay  and  data  reported  elsewhere  [34] 
generated  using  strain-specific  PCR  and  whole-genome 
restriction  enzyme  analysis. 

Rhinovirus  and  Enterovirus 

The  resequencing  results  from  HRV  and  HEV  tiles  revealed 
three  distinct  genetic  clusters  of  picomaviruses  in  the  HRV+ 
FRI  samples:  HRVA  (73),  HRVB  (five),  and  HEV  (eight; 
Table  3  and  Supplemental  Figure  2 A).  Seven  of  eight  HEVs 
were  identified  as  enterovirus  68  (EV68),  all  of  which  were 
found  in  MCRD  San  Diego  samples.  These  EV68  isolates 
shared  identical  sequences  except  for  a  single  nucleotide 
polymorphism  in  one  instance.  These  data  suggested  that 
clonal  EV68  was  endemic  in  MCRD  San  Diego  during 
2004-2005,  showing  geographic  and  temporal  clustering 
reminiscent  of  HAdV-4  outbreaks.  CAV21  was  identified  in 
one  Fort  Benning  sample,  which  corresponded  to  the 
previous  finding  in  the  Typical  FRI  sample  set,  where  two 
CAV21  were  detected  at  the  same  site.  This  might  indicate 


that  CAV21  was  transmitted  at  the  Fort  Benning  training  site. 
In  contrast,  the  78  HRV  isolates  were  identified  as  closely 
related  to  41  known  reference  serotypes  (Supplemental 
Table  3)  and  showed  broad  genetic  diversity  based  on 
phylogenetic  analysis  (Supplemental  Figure  2B).  No  domi¬ 
nant  HRV  serotype  was  identified  from  any  particular 
training  site  and  season.  This  result  suggests  HRV  coloniza¬ 
tion  in  HRV+  FRI  trainees  was  not  caused  by  endemic 
transmission  but  instead  carried  over  by  hosts.  No  geograph¬ 
ical  and/or  temporal  dependency  was  observed  for  these 
samples,  though  the  limited  sample  size  may  be  insufficient 
to  detect  significant  spatial  or  temporal  patterns  for  a  virus 
with  this  degree  of  diversity. 

Association  of  HAdV  and  HRV 

Microarray  results  showed  that  46  of  the  73  HRVA+ 
samples  from  the  HRV+  FRI  sample  set  were  also  HAdV-4+. 
The  findings  from  the  initial  sample  set  (see  “Microbial 
Distribution  in  Typical  FRI  and  Healthy  samples”)  revealed 
that  HAdV-4  may  be  negatively  associated  with  HRV.  This 
prompted  us  to  further  investigate  the  relationship  between 
these  two  vimses  in  the  HRV+  FRI  data  set.  Visually,  the 
hybridization  strengths  of  HRVA+  samples  co-infected 
with  HAdV-4  (HAdV-4+/HRVA+)  appear  relatively 
weaker  than  those  without  HAdV-4  co-infection  (HAdV-4-/ 
HRVA+;  Fig.  4).  Base  call  rates  were  compared  to  see  if 
there  was  a  significant  statistical  difference.  In  previous 
studies,  a  correlation  between  stronger  hybridization  (i.e., 
more  base  calls)  on  the  RPM-Flu  arrays  and  higher 
concentrations  of  nucleic  acid  templates  was  observed,  so 
overall  base  call  rates  indirectly  reflect  template  titers 
(unpublished  data).  The  base  call  rate  in  HRVA  tiles  for 
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HAdV4  HRVB  HRVA 


i  i  i  ii  i 


base  call  rates  -  genome  copies/jil 


HAdV4-HRV+:  33.49%*  -  5480# 


H  AdV4+HRV-K  27.88%  -  889 


HAdV4-HRV-:  1 3.0%  -  0 


*  p  <  0,007,  base  call  rate  comparison  between  HAdV4-HRV+  and  HAdV4+HRV+, 

#  p  <  0,04  viral  titer  comparison  between  HAdV4-HRV+  and  HAdV4+HRV+. 


Figure  4  HAdV-4  co-infection  effect  on  HRV  titer  in  97  HRV+  FRI 
clinical  samples.  Hybridization  profiling  in  HAdV-4  and  HRVA  tiling 
regions  shows  relative  hybridization  intensities.  HRV  titers  calculated 
with  quantitative  real-time  PCR  are  also  shown.  Average  base  call 
rates  and  genome  copies  were  calculated  based  on  three  distinct 
clusters:  HAdV-4  negative  and  HRV  positive  (HAdV-4-/HRV+), 
HAdV-4  positive  and  HRV  positive  (HAdV-4+/HRV+),  and  HAdV-4 

control  samples  and  HRVA-  clinical  samples  is  13.0%,  due 
to  background  noise.  Within  the  set  of  105  samples  chosen 
on  the  basis  of  HRV+  PCR  results,  the  average  base  call  rate 
for  HRV  tiles  in  HAdV-4-  samples  was  33.49%,  compared 
with  27.88%  for  HAdV-4+  samples  (p<0.007).  This 
suggests  that  the  titers  of  HRVA  are  lower  when  HAdV-4 
is  present. 

This  result  was  independently  validated  by  using  quanti¬ 
tative  real-time  RT-PCR.  The  average  HRV  titer  in  HAdV-4-/ 
HRVA+  samples  was  5,480  genome  copies/ pi,  approximately 
six  times  higher  than  the  880  genome  copies/ pi  seen  in 
HAdV-4+/HRVA+  samples  (/?<0.04).  This  result  may  indi¬ 
cate  that  infection  by  HAdV  interferes  with  the  replication 
efficiency  of  HRV,  or  that  HAdV  infection  is  less  likely  in 
hosts  with  higher  HRV  loads. 

Commensal  Bacteria 

H.  influenzae ,  S.  pneumoniae ,  and  N.  meningitidis  were 
detected  in  more  than  76%  of  HRV+  FRI  samples  (Table  3). 
In  addition,  sequences  generated  from  microarrays  were 
able  to  distinguish  them  at  the  strain  level.  All  85  H. 
influenzae  strains  belonged  to  nontypeable  H.  influenzae , 


negative  and  HRV  negative  (HAdV-4-/HRV-).  Analyses  of  significance 
were  generated  from  comparison  of  clusters  HAdV-4-/HRV+  and 
HAdV-4+/HRV+.  The  base  call  rate  comparison  between  HAdV-4-/ 
HRV+  and  HAdV-4+/HRV+  yielded  p <0.007.  The  viral  titer  compar¬ 
ison  of  HAdV-4-/HRV+  and  HAdV-4+/HRV+  yielded  p<0.0A.  HAdV 
human  adenovims,  HEV  human  enterovims,  HRV  human  rhino  vims 


nonencapsulated  strains  that  are  usually  harbored  as  normal 
flora  in  the  nasopharynx  [37].  The  79  S.  pneumoniae  strains 
were  classified  into  two  major  species:  Streptococcus  mitis 
and  S.  pneumoniae ,  both  of  which  belong  to  the  S.  mitis 
group  and  are  recognized  commensals  of  the  upper 
respiratory  tract  [37].  The  74  N.  meningitidis  strains  cluster 
with  members  of  serogroup  C,  which  usually  inhabits  the 
human  nasopharynx  without  causing  detectable  diseases 
[38].  The  three  commensal  bacteria  were  identified  simul¬ 
taneously  in  65%  (63/97)  of  the  HRV+  FRI  samples,  with 
no  statistically  significant  association  with  HRV  or  HAdV 
(Table  4).  Single  and  double  bacterial  colonizations  were 
more  rarely  observed.  Analysis  of  the  distributions  of 
H.  influenzae ,  S.  pneumoniae ,  and  N.  meningitidis  using 
Fisher’s  exact  statistics  showed  strongly  significant  non¬ 
randomness  (see  Table  5a;  />< 0.002).  Analysis  of  specific 
observed  and  expected  frequencies  of  the  possible  combi¬ 
nations  suggested  most  of  the  apparent  significance  derived 
from  higher-than-expected  rates  of  both  triple  infections 
(with  all  three  co-present)  and  null  infections  (with  none 
of  the  three).  Removing  the  five  null  samples  from  the 
analysis  on  the  basis  that  these  could  have  resulted  from 
poor  sample  quality  yielded  a  loss  of  significance  (see 
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Table  4  Co-colonization  of  three  major  commensal  bacterial  in  HRV+ 
FRI  samples 


HRV/HEV  positive  (84) 

HRV/HEV  negative  (13) 

Ad+  (47) 

Ad-  (37) 

Ad+  (5) 

Ad-  (8) 

Hi-Spn-Nm 

34 

23 

2 

4 

Hi-Spn 

6 

4 

0 

0 

Hi-Nm 

2 

3 

2 

0 

Spn-Nm 

1 

3 

0 

0 

Hi 

1 

2 

1 

1 

Spn 

0 

1 

0 

1 

Nm 

0 

0 

0 

0 

None 

2 

1 

0 

2 

Ad  adenovirus,  FRI  febrile  respiratory  illness,  HEV  human  enterovirus, 
HRV  human  rhino  virus,  Hi  Haemophilus  influenzae ,  Nm  Neisseria 
meningitidis ,  Spn  Streptococcus  pneumoniae 


Table  5b;  /?<0.08).  Hence,  the  data  were  suggestive  of  a 
positive  association  between  all  three.  There  is  certainly 
room  for  alternative  interpretation.  Only  five  clinical 
samples  did  not  have  any  of  the  three  commensal  bacteria. 
Three  of  them  were  infected  with  HRVA  viruses,  one  was 
infected  with  M.  pneumoniae ,  and  the  remaining  one  had 
no  identifiable  pathogens  or  commensals  at  all. 


Discussion 

Our  data  clearly  show  that  infection  with  HAdV,  the  main 
FRI  agent  in  military  training  facilities,  is  associated  with 
decreased  rates  and  titers  of  HRV,  suggesting  some  form  of 
viral  interference  between  these  groups.  The  data  also 
suggest  that  FRI-causing  pathogens  are  geographically 
segregated  and  clonal,  while  commensals  are  more  diverse 
and  do  not  display  obvious  temporal  or  spatial  clustering. 
Commensal  bacteria  are  positively  associated  with  each 
other  in  the  upper  respiratory  tract,  suggesting  some 
synergy  or  common  mechanism  favoring  colonization  in 
specific  individuals.  Furthermore,  while  the  majority  of 
picomaviruses  identified  in  this  study  were  apparently 
randomly  distributed  HRVs,  the  MCRD  San  Diego  site 
harbored  an  apparently  site-associated  strain  of  EV68  and 
Fort  Benning  harbored  CAV21,  both  of  which  were 
concentrated  among  ill  individuals,  unlike  other  identified 
picomaviruses.  Studies  should  be  undertaken  to  determine 
the  potential  health  impact  of  these  two  vimses. 

As  expected,  our  microarray  results  indicated  that  more 
than  50%  of  FRI  samples  (either  those  chosen  without 
respect  to  HRV  status  or  those  chosen  because  they  were 
HRV+)  were  infected  with  HAdV-4,  the  well-recognized 
dominant  agent  of  FRI  among  recmits.  The  data  supported 


Table  5  Statistic  analysis  of  association  of  three  commensal  bacteria  in  HRV+FRI  samples 


A.  with  5  “no  positive” 


No.  (%) 

Expected  frequency  (%) 

Observed  frequency  (%) 

Hi 

11  (11.5) 

None 

Nm 

Spn 

Spn+Nm 

None 

Nm 

Spn 

Spn+Nm 

+ 

85  (88.5) 

Hi—  0.6 

1.6 

2.2 

7.3 

Hi- 

5.2 

0.0 

2.1 

4.2 

Nm 

22  (22.9) 

Hi+  4.6 

12.1 

16.7 

56.2 

Hi+ 

5.2 

7.3 

10.4 

65.6 

+ 

74  (77.1) 

Spn 

17  (17.7) 

Expected  count  (no.):  96 

Observed  count  (No.) 

+ 

79  (82.3) 

None 

Nm 

Spn 

Spn+Nm 

None 

Nm 

Spn 

Spn+Nm 

None 

5  (5.2) 

Hi-  1 

2 

2 

7 

Hi- 

5 

0 

2 

4 

Something 

91  (94.8) 

Hi+  4 

12 

16 

54 

Hi+ 

5 

7 

10 

63 

Fisher’s  exact  test  probability  of  the  observed  distribution  (P 

-0.0021) 

B.  with  5  “no  positive” 

omitted 

No.  (%) 

Expected  frequency  (%) 

Observed  frequency  (%) 

Hi 

6  (6.6) 

None 

Nm 

Spn 

Spn+Nm 

None 

Nm 

Spn 

Spn+Nm 

+ 

85  (93.4) 

Hi-  0.3 

0.7 

1.1 

4.7 

Hi- 

0.0 

0.0 

2.1 

4.2 

Nm 

17  (18.7) 

Hi+  4.9 

10.0 

15.1 

65.9 

Hi+ 

5.2 

7.3 

10.4 

65.6 

+ 

74  (81.3) 

Spn 

12  (13.2) 

Expected  count  (No.):  91 

Observed  count  (No.) 

+ 

79  (86.8) 

None 

Nm 

Spn 

Spn+Nm 

None 

Nm 

Spn 

Spn+Nm 

None 

5  (5.2) 

Hi-  0 

1 

1 

4 

Hi- 

0 

0 

2 

4 

Something 

91  (94.8) 

Hi+  4 

9 

14 

60 

Hi+ 

5 

7 

10 

63 

Fisher’s  exact  test  probability  of  the  observed  distribution  (P 

-0.0759) 

Note.  FRI,  febrile  respiratory  illness;  HRV,  human  rhinovirus;  Hi,  Haemophilus  influenzae]  Nm,  Neisseria  meningitidis ;  Spn,  Streptococcus 
pneumoniae. 


Springer 


632 


Z.  Wang  et  al. 


previous  studies  showing  that  HAdV  outbreaks  in  training 
facilities  are  clonal  and  geographically  segregated.  The 
three  variants  identified  are  consistent  with  the  three  major 
genome  types  defined  by  restriction  enzyme  analysis  (REA) 
[34].  The  identified  geographic  pattern  of  these  three  is  also 
consistent  with  the  pattern  found  using  REA  [34], 

An  apparent  example  of  previously  unstudied  virus 
circulating  locally  in  training  sites  is  offered  by  EV68, 
one  of  the  most  rarely  reported  HEV  serotypes.  US 
surveillance  over  a  36-year  period  (1970-2005)  only 
reported  26  cases  (of  a  total  of  52,812  HEV  cases  reported), 
and  approximately  one  fourth  of  those  reports  were  in 
adults  aged  20  years  or  older  [39].  This  study  identified 
seven  EV68  cases  from  37  HRV+  FRI  trainees  exclusively 
at  MCRD  San  Diego  in  a  2-year  period  (2004-2005) 
(Table  3).  Identical  5TJTR  sequences  of  these  isolates 
revealed  by  the  microarray  agrees  with  the  monophyletic 
characteristics  of  EV68  [40].  Apparently,  the  outbreak  of 
EV68  was  geographically  segregated,  but  its  origin  and  any 
potential  disease  associations  remain  a  mystery. 

In  contrast,  multiple  HRVs  were  identified  both  in 
Healthy  and  FRI  recruits,  revealing  broad  and  scattered 
genetic  diversity  among  military  training  sites.  No  geographic 
or  temporal  clustering  of  specific  HRVs  was  observed.  Highly 
diverse  HRVs  have  also  been  documented  in  other  popula¬ 
tions  in  recent  independent  studies  [10-12].  The  clinical 
significance  of  HRV  diversity  will  be  an  interesting  and 
challenging  topic  for  future  investigation. 

Our  data  show  a  significant  negative  association  between 
HRV  and  FRI,  though  this  appears  to  be  a  secondary  effect 
of  a  negative  association  between  HRV  and  HAdV,  the 
primary  cause  of  recruit  FRI.  There  is  certainly  no  evidence 
that  HRV  is  causing  FRI  in  this  population.  Multiple  HRV 
serotypes  were  identified  at  the  same  training  site  and  even 
from  the  same  individual,  suggesting  HRV  is  ubiquitous, 
diverse,  and  commonly  subject  to  turnover. 

The  results  showed  opposing  distribution  patterns 
(negative  association)  of  HAdV  and  HRV,  with  Healthy 
samples  showing  high  rates  of  HRV  and  low  rates  of  HAdV, 
and  Typical  FRI  samples  showing  high  rates  of  HAdV  and 
low  rates  of  HRV  (Fig.  2).  This  result  suggests  some  form 
of  interference  between  HRV  and  HAdV  (or  possibly  all¬ 
cause  FRI).  This  observation  was  supported  by  results  from 
the  HRV+  FRI  samples  (Fig.  4),  which  indicated  that,  even 
among  samples  chosen  on  the  basis  of  HRV  positivity,  the 
presence  of  HAdV  is  associated  with  reduced  titers  of  HRV. 
As  all  of  these  samples  were  from  patients  with  FRI,  this 
suggests  that  the  negative  association  is  specifically 
between  HRV  and  HAdV,  rather  than  extending  to  all¬ 
cause  FRI.  Viral  interference  has  been  recently  reported  in 
two  epidemiological  studies  [41,  42],  suggesting  that  HRV 
might  have  affected  the  spread  of  pandemic  influenza  A 
(H1N1)  during  2009.  One  possible  explanation  for  this 


interference  is  that  HRV  infection  causes  the  cellular 
immune  response  to  enter  an  antiviral  state  [43]  making 
double  viral  infections  less  common  in  the  population  [42]. 
Another  possible  explanation  is  that  the  fevers  induced  by 
influenza  and  HAdV  are  high  enough  to  restrict  replication 
of  HRV. 

The  microarray  data  indicated  that  65%  of  bacteria- 
positive  samples  were  triple  colonizations,  with  fewer 
double  or  single  colonizations  occurring  than  would  be 
expected  by  chance  (Table  5).  A  recent  study  reported  that 
the  colonization  of  H.  influenzae  in  the  upper  respiratory 
tract  is  negatively  associated  with  colonization  of  either 
S.  pneumoniae  or  M.  catarrhalis,  but  the  combination  of  all 
three  shifts  the  competitive  interaction  to  positive  associa¬ 
tion  [28].  The  association  of  multiple  microbial  organisms 
may  be  more  resistant  to  clearance  by  immune  response 
[28].  Thus,  addition  of  one  more  species  might  alter  the 
competitive  balance  between  two  species  and  allow  them 
all  to  successfully  colonize  together.  Our  results  are  consistent 
with  this  finding;  however,  our  sample  size  is  too  small  to 
address  this  conclusively.  The  mechanism  of  this  multi¬ 
organism  colony  needs  to  be  further  investigated. 

From  a  clinical  standpoint,  possibly  the  most  important 
result  of  this  study  is  the  failure  to  observe  an  association 
between  colonization  with  H.  influenzae ,  S.  pneumoniae , 
N.  meningitidis ,  HRV,  or  any  combination  thereof,  with  FRI 
in  military  recruits.  The  identification  of  EV68  with  a 
specific  site,  and  the  similarity  of  the  molecular  epidemi¬ 
ology  of  this  virus  to  recognized  outbreak  pathogens  such 
as  GAS  and  HAdV,  suggests  that  this  virus  may  warrant 
further  study  as  a  potential  agent  of  significant  contagious 
disease.  The  same  is  true  for  CAV21.  From  the  standpoint 
of  microbial  ecology,  the  most  interesting  finding  was  the 
strong  negative  association,  suggestive  of  interference 
between  HRV  (whether  in  terms  of  titer  or  absolute 
occurrence)  and  HAdV-associated  FRI. 

Acknowledgments  We  acknowledge  the  efforts  of  the  entire  Naval 
Research  Laboratory  and  Naval  Health  Research  Center  team, 
especially  the  technicians  and  collection  personnel  whose  efforts  are 
represented  in  this  work.  We  also  acknowledge  the  clinic  commanders 
and  medical  staff  at  Fort  Benning,  GA;  Fort  Jackson,  SC;  Fort  Leonard 
Wood,  MO  (US  Army);  Lackland  Air  Force  Base,  San  Antonio,  TX  (US 
Air  Force);  Naval  Service  Training  Command,  Great  Lakes,  IL  (US 
Navy);  Marine  Corps  Recruit  Depot,  San  Diego,  CA;  Marine  Corps 
Recmit  Depot,  Parris  Island,  SC  (US  Marine  Corps);  and  Coast  Guard 
Training  Center,  Cape  May,  NJ  (US  Coast  Guard)  for  the  permissions, 
access,  and  assistance  necessary  to  conduct  these  studies. 

This  work  represents  NHRC  report  number  09-20,  supported  in 
part  by  the  Office  of  Naval  Research  via  Naval  Research  Laboratory 
core  funding  and  in  part  by  the  Global  Emerging  Infections 
Surveillance  and  Response  System  (GEIS),  a  Division  of  the  Armed 
Forces  Health  Surveillance  Center,  Silver  Spring,  MD  <www.afhsc.mil>, 
under  NHRC  research  work  unit  #60805.  The  views  expressed  in  this 
work  are  those  of  the  authors,  and  do  not  necessarily  reflect  the  official 
policy  or  position  of  the  Department  of  the  Navy,  Department  of 
Defense,  or  the  US  Government. 


<£)  Springer 


Microbial  Diversity  in  Sick  and  Healthy  Recruits 


633 


This  research  has  been  conducted  in  compliance  with  all  applicable 
federal  and  international  regulations  governing  the  protection  of 
human  subjects  in  research  (DoD  protocols  NHRC. 1999.0002.3 1271, 
NHRC.2003.0002,  and  NRL.2008.0013). 


References 

1 .  Russell  KL  (2006)  In:  Lenhart  MK,  Lounsbury  DE,  North  RB  Jr. 
(eds)  Textbooks  of  military  medicine:  Recruit  Medicine.  Chapter 
13,  227-253 

2.  Dudding  BA,  Top  FH  Jr,  Winter  PE,  Buescher  EL,  Lamson  TH, 
Leibovitz  A  (1973)  Acute  respiratory  disease  in  military  trainees: 
the  adenovirus  surveillance  program,  1966-1971.  Am  J  Epidemiol 
97:187-198 

3.  Van  Der  Veen  J  (1963)  The  role  of  adenoviruses  in  respiratory 
disease.  Am  Rev  Respir  Dis  88(SUPPL):  167-180 

4.  Wannamaker  LW,  Denny  FW,  Perry  WD,  Rammelkamp  CH  Jr, 
Eckhardt  GC,  Houser  HB,  Hahn  EO  (1953)  The  effect  of 
penicillin  prophylaxis  on  streptococcal  disease  rates  and  the 
carrier  state.  N  Engl  J  Med  249:1-7 

5.  Cohen  S,  Frank  E,  Doyle  WJ,  Skoner  DP,  Rabin  BS,  Gwaltney 
JM  Jr  (1998)  Types  of  stressors  that  increase  susceptibility  to  the 
common  cold  in  healthy  adults.  Health  Psychol  17:214-223 

6.  Cohen  S  (1995)  Psychological  stress  and  susceptibility  to  upper 
respiratory  infections.  Am  J  Respir  Crit  Care  Med  152:S53- 
S58 

7.  Russell  KL,  Hawksworth  AW,  Ryan  MA,  Strickler  J,  Irvine  M, 
Hansen  CJ,  Gray  GC,  Gaydos  JC  (2006)  Vaccine-preventable 
adenoviral  respiratory  illness  in  US  military  recruits,  1999-2004. 
Vaccine  24:2835-2842 

8.  Metzgar  D,  McDonough  EA,  Hansen  CJ,  Blaessing  C,  Baynes  D, 
Hawksworth  AW,  Blair  PJ,  Faix  D,  Russell  KL  Dynamic  Strain 
Turnover  and  Associated  Changes  in  Antibiotic  Resistance  Rates 
of  Group  A  Streptococcus  are  Synchronized  across  Multiple  Sites. 
In  Preparation 

9.  Metzgar  D,  Osuna  M,  Kajon  AE,  Hawksworth  AW,  Irvine  M, 
Russell  KL  (2007)  Abrupt  emergence  of  diverse  species  B 
adenovimses  at  US  military  recruit  training  centers.  J  Infect  Dis 
196:1465-1473 

10.  Lamson  D,  Renwick  N,  Kapoor  V,  Liu  Z,  Palacios  G,  Ju  J,  Dean 
A,  St  George  K,  Briese  T,  Lipkin  WI  (2006)  MassTag 
polymerase-chain-reaction  detection  of  respiratory  pathogens, 
including  a  new  rhinovirus  genotype,  that  caused  influenza-like 
illness  in  New  York  State  during  2004-2005.  J  Infect  Dis 
194:1398-1402 

11.  Kistler  A,  Avila  PC,  Rouskin  S,  Wang  D,  Ward  T,  Yagi  S,  Schnurr 
D,  Ganem  D,  DeRisi  JL,  Boushey  HA  (2007)  Pan-viral  screening 
of  respiratory  tract  infections  in  adults  with  and  without  asthma 
reveals  unexpected  human  coronavirus  and  human  rhinovirus 
diversity.  J  Infect  Dis  196:817-825 

12.  Savolainen  C,  Mulders  MN,  Hovi  T  (2002)  Phylogenetic  analysis 
of  rhinovirus  isolates  collected  during  successive  epidemic 
seasons.  Virus  Res  85:41^16 

13.  Nokso-Koivisto  J,  Kinnari  TJ,  Lindahl  P,  Hovi  T,  Pitkaranta  A 
(2002)  Human  picomavirus  and  coronavirus  RNA  in  nasopharynx 
of  children  without  concurrent  respiratory  symptoms.  J  Med  Virol 
66:417^120 

14.  Johnston  SL,  Sanderson  G,  Pattemore  PK,  Smith  S,  Bardin  PG, 
Bruce  CB,  Lambden  PR,  Tyrrell  DA,  Holgate  ST  (1993)  Use  of 
polymerase  chain  reaction  for  diagnosis  of  picomavirus  infection 
in  subjects  with  and  without  respiratory  symptoms.  J  Clin 
Microbiol  31:111-117 

15.  Wright  PF,  Deatly  AM,  Karron  RA,  Belshe  RB,  Shi  JR,  Gmber 
WC,  Zhu  Y,  Randolph  VB  (2007)  Comparison  of  results  of 


detection  of  rhinovims  by  PCR  and  viral  culture  in  human  nasal 
wash  specimens  from  subjects  with  and  without  clinical  symp¬ 
toms  of  respiratory  illness.  J  Clin  Microbiol  45:2126-2129 

16.  Hament  JM,  Kimpen  JL,  Fleer  A,  Wolfs  TF  (1999)  Respiratory 
viral  infection  predisposing  for  bacterial  disease:  a  concise  review. 
FEMS  Immunol  Med  Microbiol  26:189-195 

17.  Meyers  LA,  Levin  BR,  Richardson  AR,  Stojiljkovic  I  (2003) 
Epidemiology,  hypermutation,  within-host  evolution  and  the 
vimlence  of  Neisseria  meningitidis.  Proc  Biol  Sci  270:1667-1677 

18.  Garcia-Rodriguez  JA,  Fresnadillo  Martinez  MJ  (2002)  Dynamics 
of  nasopharyngeal  colonization  by  potential  respiratory  pathogens. 
J  Antimicrob  Chemother  50(Suppl  S2):59-73 

19.  Bogaert  D,  De  Groot  R,  Hermans  PW  (2004)  Streptococcus 
pneumoniae  colonisation:  the  key  to  pneumococcal  disease. 
Lancet  Infect  Dis  4:144-154 

20.  Hausdorff  WP,  Hajjeh  R,  Al-Mazrou  A,  Shibl  A,  Soriano-Gabarro 
M  (2007)  The  epidemiology  of  pneumococcal,  meningococcal, 
and  Haemophilus  disease  in  the  Middle  East  and  North  Africa 
(MENA)  region — current  status  and  needs.  Vaccine  25:1935- 
1944 

21.  Margolis  E,  Levin  BR  (2007)  Within-host  evolution  for  the 
invasiveness  of  commensal  bacteria:  an  experimental  study  of 
bacteremias  resulting  from  Haemophilus  influenzae  nasal  carriage. 
J  Infect  Dis  196:1068-1075 

22.  Lin  B,  Malanoski  AP,  Wang  Z,  Blaney  KM,  Long  NC,  Meador 
CE,  Metzgar  D,  Myers  CA,  Yingst  SL,  Monteville  MR,  Saad  MD, 
Schnur  JM,  Tibbetts  C,  Stenger  DA  (2009)  Universal  detection 
and  identification  of  avian  influenza  vims  by  use  of  resequencing 
microarrays.  J  Clin  Microbiol  47:988-993 

23.  Lin  B,  Blaney  KM,  Malanoski  AP,  Ligler  AG,  Schnur  JM, 
Metzgar  D,  Russell  KL,  Stenger  DA  (2007)  Using  a  resequencing 
microarray  as  a  multiple  respiratory  pathogen  detection  assay.  J 
Clin  Microbiol  45:443-452 

24.  Lin  B,  Malanoski  AP,  Wang  Z,  Blaney  KM,  Ligler  AG,  Rowley 
RK,  Hanson  EH,  von  Rosenvinge  E,  Ligler  FS,  Kusterbeck  AW, 
Metzgar  D,  Barrozo  CP,  Russell  KL,  Tibbetts  C,  Schnur  JM, 
Stenger  DA  (2007)  Application  of  broad-spectmm,  sequence- 
based  pathogen  identification  in  an  urban  population.  PLoS  ONE 
2:e419 

25.  Lin  B,  Wang  Z,  Vora  GJ,  Thornton  JA,  Schnur  JM,  Thach  DC, 
Blaney  KM,  Ligler  AG,  Malanoski  AP,  Santiago  J,  Walter  EA, 
Agan  BK,  Metzgar  D,  Seto  D,  Daum  LT,  Kruzelock  R,  Rowley 
RK,  Hanson  EH,  Tibbetts  C,  Stenger  DA  (2006)  Broad-spectmm 
respiratory  tract  pathogen  identification  using  resequencing  DNA 
microarrays.  Genome  Res  16:527-535 

26.  Wang  Z,  Daum  LT,  Vora  GJ,  Metzgar  D,  Walter  EA,  Canas  LC, 
Malanoski  AP,  Lin  B,  Stenger  DA  (2006)  Identifying  influenza 
vimses  with  resequencing  microarrays.  Emerg  Infect  Dis  12:638— 
646 

27.  Wang  Z,  Malanoski  AP,  Lin  B,  Kidd  C,  Long  NC,  Blaney  KM, 
Thach  DC,  Tibbetts  C,  Stenger  DA  (2008)  Resequencing  micro¬ 
array  probe  design  for  typing  genetically  diverse  vimses:  human 
rhino  vimses  and  entero  vimses.  BMC  Genomics  9:577 

28.  Pettigrew  MM,  Gent  JF,  Revai  K,  Patel  JA,  Chonmaitree  T  (2008) 
Microbial  interactions  during  upper  respiratory  tract  infections. 
Emerg  Infect  Dis  14:1584-1591 

29.  Santti  J,  Hyypia  T,  Halonen  P  (1997)  Comparison  of  PCR  primer 
pairs  in  the  detection  of  human  rhinovimses  in  nasopharyngeal 
aspirates.  J  Virol  Methods  66:139-147 

30.  Hendolin  PH,  Paulin  L,  Ylikoski  J  (2000)  Clinically  applicable 
multiplex  PCR  for  four  middle  ear  pathogens.  J  Clin  Microbiol 
38:125-132 

31.  Hendolin  PH,  Markkanen  A,  Ylikoski  J,  Wahlfors  JJ  (1997) 
Use  of  multiplex  PCR  for  simultaneous  detection  of  four 
bacterial  species  in  middle  ear  effusions.  J  Clin  Microbiol 
35:2854-2858 


<5)  Springer 


634 


Z.  Wang  et  al. 


32.  Malanoski  AP,  Lin  B,  Wang  Z,  Schnur  JM,  Stenger  DA  (2006) 
Automated  identification  of  multiple  micro-organisms  from 
resequencing  DNA  microarrays.  Nucleic  Acids  Res  34:5300-53 1 1 

33.  Fischer  RA  (1925)  Statistical  methods  for  research  workers. 
Oliver  and  Boyd,  Edinburgh 

34.  Kajon  AE,  Moseley  JM,  Metzgar  D,  Huong  HS,  Wadleigh  A, 
Ryan  MA,  Russell  KL  (2007)  Molecular  epidemiology  of 
adenovirus  type  4  infections  in  US  military  recruits  in  the 
postvaccination  era  (1997-2003).  J  Infect  Dis  196:67-75 

35.  Chappie  PJ  (1966)  A  survey  of  antibodies  to  Adenovirus  8  and 
Coxsackievirus  A21  in  human  sera.  Bull  World  Health  Organ 
34:243-248 

36.  Hilleman  MR,  Stallones  RA,  Gauld  RL,  Warfield  MS,  Anderson 
SA  (1957)  Vaccination  against  acute  respiratory  illness  of 
adenovims  (RI-APC-ARD)  etiology.  Am  J  Public  Health  Nations 
Health  47:841-847 

37.  Forbes  BA,  Sahm  DF,  Weissfeld  AS  (2002)  Bailey  &  Scott’s 
diagnostic  microbiology,  Eleventh  Edition.  Mosby,  New  York 

38.  Bentley  SD,  Vemikos  GS,  Snyder  LA,  Churcher  C,  Arrowsmith 
C,  Chillingworth  T,  Cronin  A,  Davis  PH,  Holroyd  NE,  Jagels  K, 
Maddison  M,  Moule  S,  Rabbinowitsch  E,  Sharp  S,  Unwin  L, 
Whitehead  S,  Quail  MA,  Achtman  M,  Barrell  B,  Saunders  NJ, 


Parkhill  J  (2007)  Meningococcal  genetic  variation  mechanisms 
viewed  through  comparative  analysis  of  serogroup  C  strain 
FAM18.  PLoS  Genet  3:e23 

39.  Khetsuriani  N,  Lamonte-Fowlkes  A,  Oberst  S,  Pallansch  MA 
(2006)  Enterovirus  surveillance — United  States,  1970-2005. 
MMWR  Surveill  Summ  55:1-20 

40.  Oberste  MS,  Maher  K,  Schnurr  D,  Flemister  MR,  Lovchik  JC, 
Peters  H,  Sessions  W,  Kirk  C,  Chatterjee  N,  Fuller  S,  Hanauer  JM, 
Pallansch  MA  (2004)  Enterovirus  68  is  associated  with  respiratory 
illness  and  shares  biological  features  with  both  the  enteroviruses 
and  the  rhinoviruses.  J  Gen  Virol  85:2577-2584 

41.  Casalegno  JS,  Bouscambert-Duchamp  M,  Morfin  F,  Lina  B, 
Escuret  V  (2009)  Rhinoviruses,  A(HlNl)v,  RVS:  the  race  for 
hivemal  pandemics,  France  2009-2010.  Euro  Surveill  14:44 

42.  Linde  A,  Rotzen-Ostlund  M,  Zweygberg-Wirgart  B,  Rubinova  S, 
Brytting  M  (2009)  Does  viral  interference  affect  spread  of 
influenza?  Euro  Surveill  14:40 

43.  Khaitov  MR,  Laza-Stanca  V,  Edwards  MR,  Walton  RP,  Rohde  G, 
Contoli  M,  Papi  A,  Stanciu  LA,  Kotenko  SV,  Johnston  SL  (2009) 
Respiratory  virus  induction  of  alpha-,  beta-  and  lambda- 
interferons  in  bronchial  epithelial  cells  and  peripheral  blood 
mononuclear  cells.  Allergy  64:375-386 


<£)  Springer 


Copyright  of  Microbial  Ecology  is  the  property  of  Springer  Science  &  Business  Media  B.V.  and  its  content  may 
not  be  copied  or  emailed  to  multiple  sites  or  posted  to  a  listserv  without  the  copyright  holder's  express  written 
permission.  However,  users  may  print,  download,  or  email  articles  for  individual  use. 


REPORT  DOCUMENTATION  PAGE 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other 
aspect  of  this  collection  of  information,  including  suggestions  for  reducing  the  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports, 
1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject 
to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  Control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO 
THE  ABOVE  ADDRESS. 


12  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

Microb  Ecol  (2010)  59^623-634 _ 

14.  ABSTRACT  (maximum  200  words) 

Military  recruits  experience  a  high  incidence  of  febrile  respiratory  illness  (FRI),  leading  to  significant  morbidity  and  lost 
training  time.  Adenoviruses,  group  A  Streptococcus  pyogenes,  and  influenza  virus  are  implicated  in  over  half  of  the  FRI 
cases  reported  at  recruit  center  clinics.  The  etiology  of  the  remaining  cases  remains  unclear.  In  this  study,  we  explore  the 
background  carriage  rates  of  rhinoviruses,  Streptococcus  pneumoniae,  Haemophilus  influenzae,  and  Neisseria 
meningitidis .  The  results  showed  that  S.  pneumoniae,  H.  influenzae,  and  N.  meningitidis  were  widely  distributed  in  86%, 
82%,  and  76%  of  recruits,  respectively.  Rhinovirus  was  identified  in  43%  of  recruits.  Asymptomatic  (healthy)  recruits  were 
as  likely  as  ill  recruits  to  be  positive  for  all  four  of  these  organisms,  suggesting  that  they  are  not  responsible  for  significant 
illness  in  training  populations.  Recruits  with  FRI  were  actually  somewhat  less  likely  to  carry  rhinovirus  than  healthy  recruits, 
the  apparent  result  of  a  negative  association  between  adenovirus  carriage  and  rhinovirus  carriage.  This  intriguing  discovery 
may  result  from  temperature  repression  of  rhinoviruses  (which  are  highly  temperature  sensitive)  by  the  fevers  induced  by 
adenovirus  infection. 


14.  SUBJECT  TERMS 

military  recruit,  febrile  respiratory  illness,  commensal,  microbial  ecolo; 

~  17.  LIMITATION  1 18.  NUMBER 

OF  ABSTRACT  OF  PAGES 

UNCL  13 

1 18b.  TELEPHONE  NUMBER  (INCLUDING  AREA  CODE) 

Standard  Form  298  (Rev.  8-98) 
Prescribed  by  ANSI  Std.  Z39-18 


18a.  NAME  OF  RESPONSIBLE  PERSON 

Commanding  Officer 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

b.ABSTRACT 

b.  THIS  PAGE 

UNCL 

UNCL 

UNCL 

1.  Report  Date  (DD  MM  YY)  2.  Report  Type 

10  04  09  Journal  Article 

4.  TITLE  AND  SUBTITLE 

Broad  Spectrum  Respiratory  Pathogen  Analysis  of  Throat  Swabs  from  Military 
Recruits  Reveals  Interference  Between  Rhinoviruses  and  Adenoviruses _ 

6.  AUTHORS 

Wang,  Z;  Malanoski,  AP;  Lin,  B;  Long,  NC;  Leski,  TA;  Blaney,  KM;  Brown,  J; 
Broderick,  M;  Stenger,  DA;  Tibbetts.  C;  Russell,  KL;  Metzgar,  D. 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Commanding  Officer 
Naval  Health  Research  Center 
140  Sylvester  Rd 
San  Diego,  CA  92106 

8.  SPONSORING/MONITORING  AGENCY  NAMES(S)  AND  ADDRESS(ES) 

Commanding  Officer  Chief,  Bureau  of  Medicine  and  Surgery 

Naval  Medical  Research  Center  7700  Arlington  Blvd 
503  Robert  Grant  Ave  Falls  Church,  VA  22042 

Silver  Spring,  MD  20910-7500 


3.  DATES  COVERED  (from  -  to) 

2000-2006 

5a.  Contract  Number: 

5b.  Grant  Number: 

5c.  Program  Element: 

5d.  Project  Number: 

5e.  Task  Number: 

5f.  Work  Unit  Number:  60805 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

09-20 

10.  Sponsor/Monitor's  Acronyms(s) 

NMRC/NMSC _ 

11.  Sponsor/Monitor's  Report  Number(s) 


